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Abstract 

Here we propose an experimental set-up in which it is possible to obtain the entan- 
T^j- | glement of a two-mode Gaussian state, be it pure or mixed, using only simple linear 

optical measurement devices. After a proper unitary manipulation of the two- mode 
Gaussian state only number and purity measurements of just one of the modes 
Q \ suffice to give us a complete and exact knowledge of the state's entanglement. 
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^ ■ 1 Introduction 

a' 



Quantum theory of information offers, in principle, a plethora of resources for 
the processing of computational and informational tasks in an efficient way, 
otherwise unattainable in the classical world [Tll2l3]H] . Initially, all the quan- 
tum communication protocols were developed for discrete systems (qubits) 
and later an equivalent formulation for the quantum continuous variable (CV) 
setting was developed. CV-systems are described by canonical conjugate op- 
erators like position and momentum or the quadrature amplitudes of the 
quantized electromagnetic field. The main reason for the development of a 
CV quantum information theory are of a practical order: the crucial steps in 
quantum communication protocols are relatively simple to implement via the 
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available experimental techniques of quantum optics [2]. Some of these tasks, 
e.g. quantum cryptography [5f6lTT] . superdense coding [8f9][T0] and teleporta- 
tion [TTlll2|ll3j . require entanglement as a key ingredient for their successful 
implementation. Therefore, a great deal of effort has been put into the qual- 
itative and quantitative study of entanglement. However, it is still a difficult 
problem to obtain in a simple experimental fashion the degree of entanglement 
of a quantum system. Practically all the useful entanglement measures devel- 
oped so far can only be computed with a complete knowledge of the quantum 
state, whose reconstruction is not a trivial experimental task. 

A direct scheme for the measurement of the entanglement of a pair of pure 
qubits, without state reconstruction, was recently proposed [H] and experi- 
mentally realized [15]. For continuous variable systems (CV-systems), being 
Gaussian states a famous example which are efficiently generated in the labo- 
ratory [2U4| . interesting experimental proposals for the measurement of entan- 
glement were presented [T6|17|I18] and implemented [19] recently. In particular, 
it was shown that by measuring the marginal (local) and global purities of a 
two- mode Gaussian state one is able to estimate the state negativity [20], al- 
though its exact value cannot be achieved via this procedure. The negativity, 
together with the entanglement of formation [5T] , are the most important en- 
tanglement measures for Gaussian states, where the latter is only analytically 
computable for symmetric Gaussian states. 

Here we present a simple experimental scheme allowing the measurement of 
the exact value of the entanglement of formation for an arbitrary symmetric 
two-mode Gaussian state. In fact, the scheme can also be adapted to obtain the 
exact value of the negativity for an arbitrary two-mode Gaussian state as well, 
being it symmetric or not. The only prior knowledge needed for entanglement 
quantification, which is also a requisite of all previous schemes, is upon the 
Gaussian nature of the two-mode state. Whether it is a pure or a mixed state 
is irrelevant for our purposes. Furthermore, even if we do not know that we 
are dealing with a Gaussian state our method can be used to experimentally 
apply the Simon separability test [22], which is a sufficient condition for the 
existence of CV-entanglement. 

A most remarkable feature of our scheme, which differentiates it from previous 
proposals, is that all the information relevant to the exact quantification of the 
entanglement of the state can be obtained, after a non-local unitary operation 
on the two-modes, via local projective measurements. After combining the 
two-modes in a beam splitter (non-local unitary operation) the entanglement 
can be quantified solely by measurements of the purity and the number of 
photons of just one of the modes. Therefore, all the information concerning 
the entanglement of the two-modes is transferred to local properties of one 
of the modes after a non-local unitary operation. It is in this sense that the 
entanglement can be seen as determined by two single types of local projective 
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measurements with no need for quantum state reconstruction or, equivalently, 
without the knowledge of the two-mode covariance matrix. 



2 Entanglement of a two-mode Gaussian state 



A two-mode Gaussian state p is completely characterized by its covariance 
matrix, with elements 7^ = (RiRj + RjRi) p — 2(Rj) p (Rj) p , where (Ri) p is 
the quantum expectation value of the observable Ri and (Ri, R 2 , R3, R4) = 
(Xi, Pi, X 2 , P2) are the quadratures of the electromagnetic field modes [2"3"j . 
Any covariance matrix can be brought via local symplectic transformations, 
i.e. without affecting its entanglement content, to the standard form [22] 



/ 



7 



C 7 T B 



where A, B, and C 7 are 2x2 real diagonal matrices given as [21] A = 
diag(n, n), B = diag(m, m), and C 7 = diag(/c x , — k p ). It can be shown [22] 
that a two-mode Gaussian state is not entangled if k x k p > 0. However, it may 
or may not be entangled when k x k p < 0. A Gaussian state is not entangled 
(separable) if, and only if [22] 

/ 1 / 2 + (l-|/ 3 |) 2 -/4>/i + / 2 , (2) 

where h = det(A), I 2 = det(B), I 3 = det(C 7 ) and I 4 = ti(AJC 7 JBJC% 'J) 
are the four invariants of the Sp(2,R) <g) Sp(2,R) group [22]. Here J is an 
anti-diagonal 2x2 matrix given by J = adiag(l, —1), det(M) stands for the 
determinant of the matrix M, tr(M) is the trace of M, and M T is the transpose 
of M. 

For a symmetric state (Ji = I 2 ) the entanglement of formation [21] can be 
written in terms of the symplectic invariants as 




E f {p) = f Jli + |/ 3 | - Jh + 2h \h\ , (3) 



where f(x) = c + (x) log 2 (c + (x)) — c_(x) log 2 (c_(x)) and c±(x) = (x" 1 / 2 ± 
x 1 / 2 ) 2 /4. There is no analytical expression of Ef for non-symmetric Gaussian 
states. However, their entanglement can be quantified by the negativity, which 
can be calculated if the above four symplectic invariants are known [To]fT7] . 
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3 The four local invariants 



It is important to note, since this is a crucial ingredient of our experimental 
proposal, that the four quantities Ji, I2, I3, and J4 are invariants by local sym- 
plectic transformations in the quadratures, or equivalently, they are invariant 
by local unitary operations in the density matrix p [22J. The knowledge of 
these four invariants allows us to completely characterize the entanglement 
properties of a two-mode Gaussian state: we can discover if it is entangled 
(Eq. (J2])) as well as how much it is entangled (Eq. (j3J)). Therefore, the mea- 
surement of these symplectic invariants is the base on which our proposal is 
built and our goal now is to present an experimental set-up in which they can 
be easily determined. 

For this purpose it is more appropriate to work with the covariance matrix V 
[25] of the creation, aj = (Xj — iPj)/y/2, and annihilation, aj = (Xj + iPj)/y/2, 
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operators for the two modes (j = 1,2) [26]. Defining v = {v\, t>2, v 3 , t> 4 ) T = 
(di, a[, a,2, cl\) T ■> t ne matrix elements of V are Vij = (—l) l+ ^{viVj + VjVi) p /2 
[27] . As we did for 7 we can represent V in terms of four block matrices of 
dimension two: 



V = t I • (4) 

\Cy V2 

In this new representation the four symplectic invariants read: J\ = det(Vi), 
J 2 = det(F 2 ), J 3 = det(CV) and J 4 = tv(V 1 ZC v ZV 2 ZC i v Z), where Z = 
diag(l, — 1). Recalling the definitions of 7^, V^, and the relation between 
Xj , Pj and ] it is straightforward to see that 

h = 4J 1; I 2 = 4J 2 , / 3 = 4J 3 , I 4 = I6J4. (5) 

Like 7, under local symplectic transformations V can be brought to the fol- 
lowing standard form [28] : 

V=\ 9 \ 6V \, (6) 



with 






Vx 



where fh s and rh c are real parameters. From now on, whenever we deal with the 
standard form of V, we will write its elements as well as any related quantity 
thereof with a tilde (~). 
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4 Experimental proposals 



We only need, then, to experimentally measure Ji, J2, J3, and J4 to completely 
quantify the entanglement of a two-mode Gaussian state. We first show a 
scheme in which we can determine the first three invariants. With these three 
invariants, as will be shown in what follows, we can obtain bounds for the 
entanglement of a two-mode Gaussian state. We then introduce two local 
unitary operations to the previous scheme, allowing us to determine J4, the 
remaining invariant. 



4-1 First scheme 

The experimental set-up necessary to measure the first three invariants is 
very simple and is depicted in Fig. [IJ It can be thought of as a simplification 
of the scheme presented in Ref. [29] to reconstruct the two-mode covariance 
matrix of a Gaussian state. Here, however, there is no state reconstruction 
and we only use linear optical devices to measure the purity and the photon 
number of the output mode a' x . Loosely speaking the linear optics apparatus 
(adjustable phase shifter and a beam splitter) can be seen as the agent re- 
sponsible for transferring the entanglement properties between the two modes 
a± and (12 to the mode a[ [26]. On the other hand, the measuring apparatus 
(photon counting and/or homodyne detection) are responsible for the local 
projective measurements which determine these properties. The phase shifter 




Fig. 1. Experimental set-up to measure three symplectic invariants (Ji, J2, and J3) 
of a two-mode Gaussian state. Input modes a± and 02 pass through a beam splitter 
(BS) with transmittance cos#. Before reaching the BS mode a\ passes through a 
phase-shifter (PS) acquiring a phase (p. After the BS measurements (M) are made 
on the output mode only. Two types of measurements are required to completely 
characterize the invariants of the input two-mode Gaussian state: the photon num- 
ber and the purity (Wigner function at the origin of the phase space) of mode 
a[. 

and the beam splitter action on the modes a x and 02 is modeled by the fol- 
lowing non-local bilinear Bogoliubov transformation [26] v' = Wv, where v' = 
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a 



\)la'(a'W and 



1) v 




(7) 

The 2x2 block matrices 1Z = diag(e 1¥ ' cos 8, e -1 ^ cos 8) and S = diag(sin 8, sin 8) 
are such that cos 8 is the transmittance T of the beam splitter and if is the 
phase shift in mode a\. The two-mode output covariance matrix is V = WVU. 
We will only need, however, the mode a[ local covariance matrix, which reads 
[26] 

v{ = n*v{R + sv 2 s* - sc ] v n - n*c v s*. (8) 

Defining J[(8,ip) = det(V{), where we explicitly write the dependence of 
det(V{) on the parameters 6 and (p, and using Eq. ([8]) we easily see that 

J! = J((0,0), J 2 = ^(|,o). (9) 

As expected Ji and J 2 are obtained when the beam splitter has, respectively, 
transmittance one (8 = 0) and reflectivity one (8 = vr/2). The determination 
of J3 is not as trivial, requiring that a set of J[ measurements be made on 
the output for distinct beam-splitter transmittances and phase-shift arrange- 
ments, as well as the measurement of the mode a[ average photon number 
N[(8,ip) = ((a')!^ + 1/2). A straightforward but tedious calculation gives 

J 3 = \(J + JV), (10) 

where 









fix 
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A^=[iV((0,0)] 2 
-2 



-1 2 



iV{(0,0)+iV{(-,0 



7T 



7T 



^(l' +JV; V4'2 



2 

7T 7T 



[ 4' 2 



11^ 



'121 



Eqs. (!9|)- (|T2|) constitute one of our central results and tell that Ji, J 2 and J 3 of 
the bipartite input state can be completely determined by measurements on 
only one of the beam-splitter output ports. This result contrasts to previous 
proposals [I~6"lll7|ll8| where measurements on both modes are always required 
to obtain the symplectic invariants, in particular J3. 

Although we do not have yet the fourth invariant, we can obtain a lower 
bound for the entanglement of formation as given by Eq. Q. Noting that J 4 
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(or equivalently I 4) is a positive quantity and that the function f(x) (Eq. (j3j)) 
is a decreasing function of x [21] we readily obtain the lower bound by setting 
the fourth invariant to zero. A similar approach allows us to derive a bound 
for the negativity [ToTfT?] . 

Remark: If, and only if, V is given as or can be brought locally to either one 
of the following two distinct forms the previous three invariants are enough to 
completely quantify the entanglement of a two-mode Gaussian state. Indeed, 
when det(CV) > one of the forms reads V\ = diag(rai, ni), V2 = diag(n2, 712), 
and Cy = diag(m s ,m*). On the other hand, if det(CV) < the matrix Cy is 
ant i- diagonal, Cy = adiag(m c , m*). Note that we cannot go from one form to 
the other via local unitary operations since we have different signs for det(CV). 
A simple calculation using these two particular covariance matrices gives 

J± = 2\J 3 \\[jJ~2- (13) 

Hence, as anticipated above, for these two cases Ji, J2, and J3 are all that is 
needed to completely characterize the entanglement of a two-mode Gaussian 
state. It is worth mentioning, however, that Eq. (fT3l) is only valid for the two 
special forms of V described above. In general, Eq. (fT3j) is no longer valid. 
Also, we have not used it in any calculations that led to the construction of 
this and the next scheme. 

4-2 Second scheme 

In order to get J4 we modify the previous scheme introducing two local unitary 
operations, as depicted in Fig. [2j The unitary transformations U\ and TJ% are 



1 -l LJ 

PS a 2 

Fig. 2. Modified experimental set-up to measure the remaining invariant J4 of a two- 
mode Gaussian state. After the local unitary transformations U\ and U2, affecting 
modes 1 and 2, respectively, the scheme is identical to the one given in Fig. [H 

chosen such that the covariance matrix describing the output modes a\ and 
a 2 is in its standard form V (See Eq. (jfJJ)). These unitary transformations 
preserve the Gaussian character of the input state and they are equivalent to 
local symplectic transformations S\ and on 7. We should bear in mind that 
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the correct transformation is dependent on the one-mode covariance matrices 
Vi and V 2 , which fortunately can be determined locally by standard homodyne 
detection. Moreover, the important point here is that S\ and S 2 (or U\ and 
U 2 ) always exist and that they are local symplectic transformations. In the 
Appendix we show how to express these transformations as a function of the 
input local covariance matrices V\ and V 2 . 

In the standard form V a direct calculation shows that the fourth invariant is 
given as 

J 4 = 2nin 2 (\fh s \ 2 + |m c | 2 ) . (14) 

Our task then reduces to the determination of hi, h 2 , \hi s \ 2 , and \m c \ 2 mea- 
suring only the output mode a[, the one obtained after d\ and d 2 enter the 
beam splitter. 

We first note that h\ and h 2 are trivially related to invariants already measured 
in the previous scheme: h\ = y/Ji and h 2 = y/J^- Second, J 3 = |m s | 2 — \m c \ 2 
is also known, implying that we just need to measure either |m s | 2 or \fh c \ 2 to 
obtain J4. Analyzing V\ we see that, 



\m r 



Re (rh. 



Im (m, 



hi + h 2 
2 

hi + h 2 



mo), 


(15) 


TV \ 


(16) 


4' ) ' 


7T 7T\ 


(17) 


4' 2J ' 



where J[(9,ip) = det(V{) and N[(9,ip) = ((a')ia[ + 1/2) is mode a[ average 
photon number, all quantities determined for a given set of parameters 9 and 
if. Note that although in the standard form the matrix elements m s and 
m c are real, we have assumed them to be complex, which simplifies the unitary 
operations U\ and U 2 necessary to transform V to V. See the Appendix for 
more details. 



Looking at Eqs. (|T5|) - (|T7|) we see that, similar to the first scheme, only two 
types of measurements are needed to determine J 4 , namely J[(9, cp) and N[(9, (p). 

Additionally, the scheme just presented can be employed without relying on 
the first one. Indeed, assuming that we are always dealing with an experimental 
set-up as depicted in Fig. [21 the first two invariants are 



J2 







1 2 



(18) 
(19) 
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With the aid of Eqs. (fT5]) - (fTT|) . the last two invariants, J3 and J4, are readily 
obtained. 

In summary, given a general two-mode Gaussian state the first scheme allows 
us to determine the first three invariants (J\, J2, and J 3 ) while the second 
scheme gives us the fourth as well as the other three invariants. With these 
four invariants the entanglement content of a two-mode Gaussian state is fully 
determined. 



5 Experimental feasibility 

Our last task is to explain how J[ and N[ (or equivalently J[ and N[ ) are 
obtained from experimentally measurable quantities. Firstly, N[ is simply the 
output mode a[ photon number (or intensity) and can be easily determined by 
a photodetection process. The determinant J[, on the other hand, is connected 
to the purity of mode a[ by the following expression [ToTTl?] : Tr{(p[) 2 } = 
l/(2y[j[). Moreover, one can prove that [2] W(0) = 1/{2ti^j[) where W(0) 
is the Wigner function of mode a[ at the origin of the phase space. Therefore, 
any technique developed to measure the purity and/or W(0) can be employed 
to determine J[. Two interesting proposals were presented in Refs. [30][T8] . 
The first one [20], implemented in Ref. [21], shows that the photon counting 
statistics allows one to obtain W(0) without any sophisticated data processing. 
The second one [18] , implemented in Ref. [32], employs only a tunable beam 
splitter and a single-photon detector to obtain the purity of mode a' v Remark 
that both schemes do not require homodyne detection and permit a direct 
access to N[ as well. However, homodyne detection [33] can also be employed 
for each set of the parameters 6 and <p to reconstruct the output mode a[ 
covariance matrix, and thus leading to the immediate calculation of J[ and 
N[. Nevertheless, a complete single-mode reconstruction is more than we need 
to fully characterize the entanglement of a two-mode Gaussian state and also 
more experimentally demanding than the previous two techniques. 

In fact what determines the detection scheme to be employed is the spe- 
cific scheme detection efficiency in the light frequency range and intensity in 
use. Obviously a photon counting scheme for the measurement of the Wigner 
function at the origin of the phase space is more appropriate for low inten- 
sity light fields. On the other hand, homodyne reconstruction is more ap- 
propriate for continuous bright light fields, unless an alternative procedure 
allowing access to the photon statistics is implemented, such as for pulsed 
fields detected by on/off avalanche photodetectors [M]. Photon counting (im- 
plemented via avalanche devices with saturated gain where each photon pro- 
duces a detectable signal (current) at the output) can also be simulated by 
dual-homodyne schemes [33] . Such schemes may be useful in infrared and op- 
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tical frequencies since the present day avalanche photodiodes possess very low 
efficiencies at such communication frequencies [35] . 

Whichever the scheme employed imperfect detection and signal losses blur the 
measurement outcomes avoiding an exact determination of the invariants as 
delineated above. However, in any situation here considered, the losses can 
be modeled by a beam splitter of transmittance t] followed by an ideal de- 
tector, being 77 attributed to the overall efficiency of the detection scheme 
[32] . Thus, the actual measured quantities can be corrected by an appropriate 
rescaling. Firstly, for homodyne detection the procedure to measure the de- 
tection efficiency r\ hom is well established from squeezing experiments [3U|I3"2"] . 
Since J[ = I'JA = (V{ tmin V{ jmax ) /4 and N[ = (y{ >min + V{ m J) /2, where 
Vlmin an d V{max are the squeezed and anti-squeezed quadratures variances, 
respectively, one can show that in terms of the actual measured variances 
Kmm and V lmax they are corrected to V{ jmin = (v( )Tni „ - 1 + r) hom ) /r) hom 
and V{ max = (yi max - 1 + Vhom) /Vhom [22] • Secondly, the determination of 
J[ and N[ via photocounting techniques relies on single photon detection of 
output mode 1 (given by a[) after the two input modes are recombined in a 
beam splitter of transmittance T = cos 9. The detector efficiency 77 is modeled 
by a beam splitter of transmittance 77 followed by an ideal detector, which in 
Refs. [T8H32] responds with only two measurement outcomes, while in Refs. 
[30|3T] is a number resolving photodetector. In either case the effects of non- 
ideal detection can be simply taken into account by substituting T — > 77T, 
whose overall effect is to rescale the proper transmittance necessary to obtain 
all the invariants 



6 Conclusion 

We have shown that continuous variable entanglement can be directly detected 
and quantified with few simple measurements. In particular, we have proposed 
an experimental scheme in which only two types of measurements, i.e. purity 
(Wigner function at the origin of the phase space) and photon number of 
a single-mode, are required to characterize the entanglement of a two-mode 
Gaussian state. Our proposal is valid for either pure or mixed states as well 
as for either symmetric or non-symmetric Gaussian states, giving always the 
exact value of the state's amount of entanglement. Furthermore, our scheme 
can also be seen as a simple procedure to obtain all the four independent 
symplectic invariants of a covariance matrix. This allows one to easily test 
for the existence of entanglement even in non-Gaussian states via the Simon's 
sufficient condition of inseparability. Finally we remark that the same scheme 
could be appropriately adapted for analyzing continuous variable entangle- 
ment in other systems than light fields, such as for the motional degrees of 



10 



freedom of trapped ions [37] , or even for Bose-Einstein condensates [38] . This is 
the case whenever similar processes to the beam splitter and the phase shifter 
operations can be applied and whenever proper local projective measurements 
or local state reconstruction can be implemented. 



Note: After completion of this manuscript we devised a fully local procedure 
by which we can measure the four invariants of a general two-mode Gaussian 
state [39] • Contrary to the schemes here presented, in Ref . [39] there is no need 
for a non-local unitary operation (as given here by the beam splitter). The 
trade-off is the inclusion of a classical communication channel and one-mode 
parity measurements. 
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A Obtaining the standard form V 



The covariance matrix V is generally written as: 



V 



Vi C v 

c v v 2 



' Tli m l m s m c ^ 

ml rii m* c m* 
m* m c ri2 m2 
ym* m s m* 2 n 2 J 



(A.1) 



where Vj, j = 1, 2, and Cy are block matrices of dimension two. Unitary local 
operations preserving the Gaussian character of a state described by V are 
mapped to the following local symplectic transformation 



S 




(A.2) 
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where Sj, j = 1, 2, is given as 



/ g lOtj C0 g] 1 ^ gift gjj^ ^. 



e sinh 9j e lolj cosh 



(A.3) 



with aij, /3j, and 9j real parameters. The first one is related to rotations of the 
quadratures of the quantized electromagnetic field and the last one to local 
squeezing operations. The new covariance matrix V is connected to V by the 
following relation [40PTT] . 

V = SVS\ (A.4) 

which implies that 



1} Sj\)S]. (A.5) 

C v = SiCySt (A.6) 

Explicitly, Eq. (|A.5[) gives for the non-diagonal term of Vf. 

rhj = e~ 2iaj mj cosh 2 9j + e 2i/3j m* sinh 2 9j 

+e i(/3,-^) n . sinh (2^.). (A.7) 

By setting rrij = |mj|e w and solving for rhj = we get 



aj = ft = *±1., (A.8) 
tanh(20 7 ) = J-^i. (A.9) 

Note that this is one of several other solutions. However, for our purposes, one 
is enough to prove that V can be locally transformed to V. 
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